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Abstract 

Background: Little is known about risk factors for severe outcomes in patients infected with H5N1 and no systematic review 
has been conducted. Understanding risk factors is an important step for prioritizing prophylaxis or treatment in the event of 
a pandemic. 

Objectives: lo systematically evaluate risk factors for severe outcomes in patients with avian influenza H5N1 infection. 
Data sources: MEDLINE, EMBASE, CINAHL, GlobalHealth, and CENTRAL through March 2011 

Eligibility criteria for selecting studies: Observational studies of any design published in English, French, Spanish, German 
or Korean that reported on risk factor-outcome combinations of interest in participants with confirmed H5N1 infections. 
Outcomes considered included death, ventilator support, hospital and ICU admission, pneumonia, and composite 
outcomes. 

Study appraisal: Risk of bias was assessed using the Newcastle-Ottawa scale (NOS). 

Results:\Ne identified 20 studies reporting on 999 patients infected with H5N1 . The majority of studies (n = 14, 70%) were at 
intermediate risk of bias, i.e. 4-6 points on the NOS. Females were at increased risk of death (OR 1 .75, 95% CI 1 .27-2.44), 
while young age, in particular <5 years of age (OR 0.44, 95% CI 0.25-0.79 for death), was protective. Data on traditional risk 
factors was scarce and requires further studies. Another major limitation in the published literature was lack of adjustment 
for confounders. 

Interpretation: Females were at increased risk for complications following H5N1 infection while young age protected 
against severe outcomes. Research on traditional risk factors was limited and is required. 
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Introduction 

Outbreaks of highly pathogenic avian influenza H5N1 infec- 
tions in poultry were first reported from Guangdong, China in 
1996 [1]. Since a cluster in Hong Kong in 1997, transmission to 
humans has been intermittently reported and, as of April 2013, a 
total of 628 confirmed cases with 374 deaths of H5N1 infected 



patients (59.6% case fatality rate) have been reported to the World 
Health Organization (WHO) from 15 countries [2]. 

Understanding risk factors for severe outcomes is an important 
step in order to prioritize prophylaxis or treatment in the event of 
a pandemic. As part of a systematic review sponsored by the 
WHO, we aimed to identify risk factors for severe outcomes or 
complications of influenza infections. We report our findings 
specific to highly pathogenic avian influenza H5N1. 
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Records identified through 
database searching: 
EMBASE 
MEDLINE 
WORLD HEALTH 
CINAHL 
CENTRAL 

z 



n = 
n = 
n = 
n = 
n = 



30,501 
20,579 
7,353 
4,502 
602 



Additional records identified 
through other sources 
From bibliographies of relevant 
studies and review articles; n = 16 



- of records screened after duplicates removed 
n = 47,874 



- of records meeting eligibility criteria for full text 
screening; n = 1,817 



# of records excluded based on screening of 
titles and abstracts; n - 46,057 



Studies excluded in full-text screening/data collection: 






No original data provided 


n = 


226 


No comparison of interest reported 


n = 


732 


Not meeting diagnostic criteria 


n = 


193 


Other reasons for exclusion 


n = 


25 


Language other than English, French, German, Spanish 


n = 


270 


# of articles excluded 


n = 


1446 


Ecological design 


n = 


117 


# Total articles excluded 


n = 


1563 





# of articles with data collected 
n = 254 



# of articles on H5N1 
Avian Influenza 
n = 20 



- of articles on non-H5Nl 
n = 234 
(not reported here) 



Figure 1. Flow diagram of included and excluded studies. 

doi:1 0.1 371 /journal.pone.0089697.g001 



Materials and Methods 

The methods were reported in detail previously [3] . Briefly, we 
included studies reporting on at least one risk factor-outcome 
combination of interest irrespective of the study design. Only 
single-case reports were excluded, but case series reporting on 
patients with and without a specific outcome for a particular risk 
factor of interest were included. Articles based on secondary 
analyses of previously published data and review articles were not 
eligible. Outcomes of interest included community-acquired 
pneumonia, mortality, hospitalization, intensive care unit (ICU) 
admission, need for ventilator support, and composite outcomes. 
Studies in English, French, German, Spanish, and Korean were 
included. 

We searched MEDLINE, EMBASE, CINAHL, Global Health, 
and the Cochrane Central Register of Controlled Trials (CEN- 
TRAL) up to March 29 of 2011 and reference lists of identified 
articles and review articles. Screening and data abstraction were 
conducted independently and in duplicate by pairs of reviewers 
using piloted and standardized forms. Two reviewers assessed the 
risk of bias independently using the Newcasde-Ottawa scale (NOS) 
[4]. 



A random-effects model was used to calculate summary 
estimates using Review Manager 5.0 (Cochrane Collaboration) 
[5]. Risk estimates are reported as odds ratio (OR) with 95% 
confidence intervals (CI). If there was a potential overlap in study 
populations across studies in a meta-analysis, studies were 
excluded in a sensitivity analysis in order to minimize any effect 
of such overlap. The I 2 statistics was used to evaluate heterogeneity 
[6] . Two a priori defined subgroup analysis based on the place of 
enrollment (community versus hospital versus ICU population) 
and risk of bias in included studies were conducted in the presence 
of significant heterogeneity defined by I 2 >60% [7]. Funnel plots 
were visually interpreted to assess publication bias but no formal 
statistical tests were conducted due to the small number of studies 
for the vast majority of risk factor-outcome combinations. 

The World Health Organization (WHO) funded the study. The 
protocol was reviewed by the WHO and the WHO's suggestion 
were incorporated into the protocol. 

Results 

Out of 47,874 tides and abstracts screened, 20 studies 
[8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27] that 
included 999 patients reported risk factors for severe outcomes for 
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male sex 


female sex 




Odds Ratio (Non-event) 


Study or Subgroup Events 


Total 


Events 


Total 


Weight 


M-H, Random, 95% CI 


0.1.1 Community 












Fasina 2010 3 


32 


24 


53 


6.3% 


8.00 [2.17, 29.53] 


Kandun 2008 47 


63 


56 


64 


12.4% 


2.38 [0.94, 6.06] 


Olsen 2005 6 


10 


10 


13 


3.3% 


2.22 [0.36, 13.54] 


Pollack 1998 2 


8 


4 


10 


2.6% 


2.00 [0.26, 15.38] 


Tran 2004 5 


6 


3 


4 


1.1% 


0.60 [0.03, 13.58] 


Subtotal (95% CI) 


119 




144 


25.7% 


2.95 [1 .54, 5.63] 


Total events 63 




97 








Heterogeneity: Tau 2 = 0.00; Chi 2 = 3.71 , df = 


4 (P = 0.45); I 2 = 


0% 




Test for overall effect: Z = 3.27 (P = 0.001 ) 










0.1.2 Hospitalized patients 












Chotpitayasunondh 2005 6 


8 


2 


4 


1.7% 


0.33 [0.03,4.19] 


Giriputro 2008 10 


13 


11 


14 


3.3% 


1.10 [0.18, 6.76] 


Hien 2004 5 


6 


3 


4 


1.1% 


0.60 [0.03, 13.58] 


Kandun 2006 1 


5 


3 


3 


0.9% 


21 .00 [0.64, 689.99] 


Liem 2009 14 


37 


12 


30 


11.0% 


1.10 [0.41, 2.94] 


Oner 2006 1 


3 


3 


5 


1.2% 


3.00 [0.15, 59.89] 


Peiris 2004 1 


2 


0 


0 




Not estimable 


Sedyaninhsih 2007 22 


33 


19 


21 


4.1% 


4.75(0.93,24.17] 


Soepandi 2010 7 


8 


11 


14 


1.8% 


0.52 [0.05, 6.09] 


WHO 2007 71 


129 


83 


127 


42.4% 


1.54 [0.93, 2.55] 


Yu 2007 2 


3 


3 


3 


0.8% 


4.20 [0.12, 151.97] 


Yu 2008 7 


11 


10 


15 


4.0% 


1.14 [0.22, 5.84] 


Yuen 1998 2 


5 


3 


7 


2.0% 


1.13 [0.11, 11.60] 


Subtotal (95% CI) 


263 




247 


74.3% 


1.47 [1.00, 2.14] 


Total events 149 




163 








Heterogeneity: Tau 2 = 0.00; Chi 2 = 7.72, df = 


11 (P = 0.74); I 2 


= 0% 




Test for overall effect: Z = 1 .97 (P = 0.05) 










Total (95% CI) 


382 




391 


100.0% 


1.75 [1.26, 2.43] 


Total events 212 




260 








Heterogeneity: Tau 2 = 0.00; Chi 2 = 14.79, df 


= 16(P = 


0.54); I 


2 = 0% 




Test for overall effect: Z = 3.35 (P = 0.0008) 










Test for subgroup differences: Chi 2 = 


3.32, df = 1 (P = 


0.07), I 


2 = 69.8% 



Odds Ratio (Non-event) 
M-H, Random, 95% CI 



0.001 



0.1 1 10 1000 
male sex female sex 



Figure 2. Forest plot comparing mortality in females as compared to males. 

doi:1 0.1 371 /journal.pone.0089697.g002 

Table 1. Summary estimates of risk-factor-outcome comparisons during H5N1 Avian influenza with odds ratio (95% confidence 
interval), I 2 , and number of included studies (n); where applicable, sensitivity analysis in italic fonts. 







Pneumonia 


ICU admission 


Ventilator support 




All-cause mortality 


Female sex 


3.03 (0.91-10.00), 0%, n=4 


3.70 (0.33-50.00), n/a, n = 1 


2.94 (0.43-20.00), 0%, n 


= 3 


1 .75*0 .27-2.44), 0%, n = 1 7 




5.00 (0.91-33.33), 0%, n = 2 








2.27*11.12-4.55), 50%, n = 4 


Elderly vs. non elderly adults 


n/a 


n/a 


n/a 




1.04 (0.12-9.08), 0%, n = 2 


Paediatric vs. adults 


0.16*(0.03-0.86), 0%, n = 3 


0.50 (0.02-12.90), n/a, n = 1 


0.62 (0.09-4.45), 0%, n = 


3 


0.72 (0.34-1.51), 71%, n = 14 




0.12 (0.01-1.03), 0%, n=2 








0.21 (0.03-1.41), 87%, n = 4 


Paediatric vs. 


0.12 (0.01-1.03), 0%, n = 2 


0.50 (0.02-12.90), n/a, n = 1 


0.93 (0.10-8.27), 0%, n = 


2 


0.48 (0.13-1.79), 73%, n = 6 


non-elderly adults 










0.25 (0.04-1.55), 87%, n = 3 


<5 y in paediatric 


0.27 (0.01-7.59), 50%, n = 2 


n/a 


0.1 1 (0.00-3.35), n/a, n = 


1 


0.44*(0.25-0.79), 0%, n = 7 


population 


0.06 (0.00-1. 23),n/a, n=1 








0.42*(0. 18-0.97), 0%, n=2 


<2 y in paediatric population 


0.71 (0.02-22.34), n/a, n = 1 


n/a 


n/a 




1.13 (0.03-37.44), n/a, n = 1 


Any risk factor or 


1.71 (0.35-8.36), 0%, n = 2 


0.60 (0.05-6.79), n/a, n = 1 


n/a 




1.04 (0.20-5.35), 9%, n = 2 


co-morbidity 


1.43 (0.18-1 1.1 ),n/a, n = 1 










Smoking 


9.00 (0.30-271.65), n/a, n = 1 


5.00 (0.15-166.59), n/a, n = 1 


0.58 (0.06-5.55), 0%, n = 


2 


0.56 (0.05-5.55), 0%, n = 2 



*: statistical significance, ICU: admission to intensive care unit, y: years; sensitivity analysis in italic: exclusion of studies with a potential overlap in patient population. 
doi:1 0.1 371 /journal.pone.0089697.t001 
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age <5 years 5-1 8 years 



Odds Ratio 



Studv or SubarouD 


Events 


Total Events 


Total 


Weight 


M-H, Random, 95% C 


6.1.1 Community 














Olsen 2005 


1 


2 


5 


5 


2.5% 


0.09 [0.00, 3.59] 


Pollack 1998 


1 


7 


1 


4 


3.5% 


0.50 [0.02, 11.09] 


WHO 2006 


9 


21 


49 


Ol 


oc no/ 




Subtotal (95% CI) 




30 




90 


42.1% 


0 44 TO 18 1 091 


Total events 


11 




55 








Heterogeneity: Tau 2 = 


0.00; Chi 2 = 


0.76 df = 


o CP = 


KJ.OO), 


|2 - no/ 
I — U /o 




Test for overall effect: Z = 1 .77 (P 


= 0.08) 










6.1.2 Hospitalized patients 












Oner 2006 


0 


2 


4 


ft 


9 Q% 

£.. \3 /O 


n 1 1 rn nn ^ ^si 

U. I 1 [u.UU, O.tJOJ 


Sedyaninhsih 2007 


4 


5 


17 


24 


6.1% 


1.65 [0.16, 17.47] 


WHO 2007 


12 


27 


69 


105 


46.1% 


0.42(0.18,0.99] 


Yuen 1998 


1 


5 


1 


z 


O TO/ 


u.zd [u.un, o.ooj 


Subtotal (95% CI) 




39 




137 


57.9% 


0 44 TO 20 0 951 


Total events 


17 




91 








Heterogeneity: Tau 2 = 


0.00; Chi 2 = 


1.94, df = 


3(P = 


058); 


I 2 = 0% 




Test for overall effect: Z = 2.10 (P 


= 0.04) 










Total (95% CI) 




69 




227 


100.0% 


0.44 [0.25, 0.79] 


Total events 


28 




146 








Heterogeneity: Tau 2 = 


0.00; Chi 2 = 


2.70, df = 


6(P = 


0.85); 


I 2 = 0% 





Odds Ratio 
M-H, Random, 95% CI 



Test for overall effect: Z = 2.74 (P = 0.006) 

Test for subgroup differences: Chi 2 = 0.00, df = 1 (P : 



0.99), I 2 = 0% 



0.001 0.1 1 10 1000 
5-1 8 years age <5 years 



Figure 3. Forest plot comparing mortality in children less than 5 years of age to 5-18 years of age. 

doi:1 0.1 371 /journal.pone.0089697.g003 



H5N1 avian influenza (Figure 1). All studies were cohort studies 
and published in English between 1998 and 2010. Evidence of 
H5N1 influenza infection was by laboratory confirmation in all 
studies. Only two studies reported on cases outside of Asia. 

As assessed by the Newcasde-Ottawa scale, the majority of 
studies were at intermediate risk of bias (i.e. 4-6 points; n= 14, 
70%) while 4 (20%) were at highest risk of bias (i.e. 3 or less points) 
and 2 (10%) at lowest risk of bias (i.e. 7-9 points). Although the 
adjusted odds ratios were reported in only one study, the use of 
adjusted odds ratios instead of the crude odds ratios would not 
have changed our assessment of the risk factors. 

Female sex was associated with higher all-cause mortality (OR 
1.75, 95% 1.27-2.44) (Figure 2, Table 1). We also found a non- 
significant trend towards higher rates of pneumonia (OR 3.03, 
95% CI 0.91-10.00), ICU admission (OR 3.70, 95% CI 0.33- 
50.00), and need for ventilator support (OR 2.94, 95% CI 0.43- 
20.00) in women. Because overlap in study populations was likely 
across some studies for mortality and pneumonia outcomes, a 
sensitivity analysis was conducted excluding studies with a 
potential overlap. The findings were similar (Table 1). 

We found no evidence that older age was a risk factor for more 
severe outcomes when comparing >65 years of age to non-elderly 
adults (OR for all-cause mortality 1.04, 95% CI 0.12-9.08). 
Children < 1 8 years of age were at lower risk of developing 
pneumonia (OR 0.16, 95% CI 0.03-0.86) when compared to 
adult patients. There was a non-significant tendency towards lower 
mortality rates (0.48, 95% CI 0.13-1.79) and lower risk of 
pneumonia (OR 0.12, 95% CI 0.01-1.03) when compared to non- 
elderly adults. We found that children of <5 years of age were at 
lower risk of death when compared to older children 5-18 years of 
age (OR 0.44, 0.25-0.79; Figure 3) with a non-significant trend 



towards a lower likelihood of developing pneumonia or need 
ventilator support. 

Only two studies reported on the presence of co-morbidities as 
potential traditional risk factors. Although no significant differ- 
ences were found, there was a non-significant tendency towards 
increased risk of pneumonia in the presence of co-morbidities (OR 
1.71, 95% CI 0.35-8.36). No significant differences in outcomes 
were found for smoking. 

The only comparison with significant heterogeneity was for all- 
cause mortality when comparing paediatric patients to non-elderly 
adults (I 2 = 73%). A subgroup analysis excluding studies at highest 
risk of bias (Newcastle-Ottawa scale of 3 or lower) reduced the I 2 
to 0% with a smaUer overaU effect (OR 0.90, 95% CI 0.51-1.57 
versus OR 0.48, 95% CI 0.13-1.79) in the primary analysis. 

Discussion 

We found that female sex was associated with worse outcomes 
following H5N1 infection and that children, particularly those <5 
years old, had better outcomes. Data on underlying co-morbidities 
as risk factors for severe outcomes were scarce and therefore there 
is no evidence demonstrating that co-morbidities affect the 
outcome in patients infected with H5N1. 

At the time of our literature search, 539 human cases with 
H5N1 infections from 15 countries had been reported to WHO 
[28] . In the 20 studies identified in our systematic review, a total of 
999 cases were reported. While it seems unlikely that a large 
number of confirmed H5N 1 infections have not been reported to 
WHO [29], the excess in cases in combination with identical time 
periods and regions reported in included studies indicates overlap 
in patient populations reported in the literature. In particular, 
there was an overlap between two large WHO reports (Appendix 
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A, references 1 1 and 1 3) and studies published by local 
investigators. In a sensitivity analysis we excluded studies that 
appeared to have included some participants reported in another 
included study, the results were in keeping with the primary 
analyses. Since the overlap in participants between studies was 
never 100%, the full dataset was presented as the primary analysis. 

Female sex was associated with a higher mortality risk as well as 
with a non-significant trend for other severe outcomes. Our 
findings are in line with a report from Egypt where the H5N1 case 
fatality rate was 49% in females compared to 8% in males [30]. Of 
note, the overall mortality rate was lower in Egypt as compared to 
mortality rates reported from other countries (34 versus 66%) [30]. 
A higher mortality in females was also found in a review of 294 
cases occurring between 2006 and 2010 published by Fiebig et al. 
[29]. Women infected with H5N1 were older [29] and most of the 
affected women were housewives taking care of poultry [31]. In 
contrast, the males were largely pre-school boys that got exposed 
while playing outdoors [31]. An exposure to a greater inoculum is 
one possibility for the worse outcomes in women. Another is that 
the sex differences were confounded by age: females were likely 
exposed at an older age than boys, and it was shown that younger 
age, in particular <5 years of age, was associated with a lower risk 
for severe outcomes. This is in keeping with the findings in the 
above mentioned study by Fiebig et al. There, the authors found 
an even more pronounced difference with five to six times lower 
odds for death for children up to 9 years of age as compared to a 
number of other age groups [29]. 

In contrast to non-H5Nl influenza infections where older age 
was shown to be a risk factor for severe outcomes for seasonal 
influenza as well as during the 2009 H1N1 pandemic [3], this does 
not seem to be true for H5N1 infections. One can hypothesize that 
previous exposure to other subtypes may have offered protection 
against H5N1 [32]. 

One limitations of this systematic review is the potential for 
reporting bias. In particular asymptomatic patients or patients 
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with minor symptoms were likely underrepresented in the cases 
reported. In fact, a 1 to 2% infection rate in those exposed was 
found in serological studies which would result in much higher 
number of cases than reported [33], which in turn is regarded to 
be an overestimate by others [34]. Adjusted risk estimates also 
were not available. Thus, modifying risk factors such as differing 
time to oseltamivir treatment or the availability of health care may 
have biased outcomes [29,35]. Also, we did not consider articles in 
Chinese or other languages spoken in the regions where H5N1 is 
endemic. However, we did not find any relevant articles on H5N1 
when screening titles and abstracts of the 50 Chinese articles that 
were excluded. Among the remaining 220 articles that were 
excluded due to language, we only identified three review articles 
on H5N1 when searching for the key words "avian" or "H5N1" in 
the tide. Thus, we believe that we did not miss any relevant articles 
on Avian Influenza in the databases that we searched due to 
restrictions on language. Finally, sample sizes for some rarely 
reported risk factor-outcome comparisons were small explaining 
the low precision and likely the lack of statistical significance for 
some of these trends that we noted. 

In summary, we found evidence that females were at increased 
risk for complications following H5N1 infection while young age 
protected against severe outcomes. Research on traditional risk 
factors is scarce and is required. 
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